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C\( ' ABSTRACT 

Aims. We attempt to study spatial variations of the spectrum and rotation measures (RMs) of the large diameter, high-latitude super- 
nova remnant (SNR) CTA 1 . 

Methods. We conducted new A6 cm and All cm observations of CTA 1 using the Urumqi 25-m and Effelsberg 100-m telescopes. 
Data at other wavelengths were included to investigate the spectrum and polarisation properties. 

Results. We obtained new total intensity and polarisation maps at A6 cm and Al 1 cm with angular resolutions of 9'.5 and 4'4, respec- 
tively. We derived a spectral index of a = -0.63 + 0.05 (S v oc y a ) based on the integrated flux densities at 408 MHz, 1420 MHz, 
2639 MHz, and 4800 MHz. The spectral index map calculated from data at the four frequencies shows a clear steepening of the 
spectrum from the strong shell emission towards the north-western breakout region with weak diffuse emission. The decrease of the 
spectral index is up to about Aa = 0.3. The RM map derived from polarisation data at A6 cm and All cm shows a sharp transition 
between positive RMs in the north-eastern and negative RMs in the south-western part of the SNR. We note a corresponding RM 
pattern of extragalactic sources and propose the existence of a large-diameter Faraday screen in front of CTA 1, which covers the 
north-eastern part of the SNR. The RM of the Faraday screen is estimated to be about +45 rad m -2 . A RM structure function of CTA 1 
indicates a very regular magnetic field within the Faraday screen, which is larger than about 2.7 fiG in case of 500 pc distance. 
Conclusions. CTA 1 is a large-diameter shell-type SNR located out of the Galactic plane, which makes it an ideal object to study 
its properties without suffering confusion. The rare breakout phenomenon known for CTA 1 is confirmed. We identify a Faraday 
screen partly covering CTA 1 with a regular magnetic field in the opposite direction to the interstellar magnetic field. The detection 
of Faraday screens in the Galactic plane is quite common, but is difficult at high latitudes where the polarisation angles of weak 
polarised background emission are rotated. RMs from extragalactic sources are needed for this purpose, although the number density 
of extragalactic RMs is still small despite of significant observational progress. 
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1- Introduction the low-density region towards the north-west, where the mag- 
netic field is not sufficiently compressed and the limb-brightened 
The supernova remnant (SNR) CTA 1 (Gl 19.5+10.2) was dis- shell is not formecL Breakout phenomena of SNRs can be used 
O ; covered by Harris & Roberts (1960). Since then extensive ob- t0 study the non -uniformity of the interstellar medium. Spatial 
servations at radio, optical, and X-ray bands have been made, variations in the spectral index could be a diagnosis of the par- 
. Radio maps at high-angular resolution were observed with the ticle acce leration mechanisms. The steepening of the spectrum 
Effelsberg 100-m telescope at 1720 MHz and 2695 MHz (Sieber towards the diffuse region was qualitatively shown by Pineault 
, et al. 1979, 1981) and with the DRAO synthesis array combined et al (1997) based on 408 MHz and 1420 MHz data. However, 
with data from the Effelsberg 100-m telescope at 408 MHz and observations at higher frequencies are needed to obtain a more 
1420 MHz (Pineault et al. 1993, 1997). CTA 1 shows a well- precise spatial d i stribut ion of the spectral index. 
Fz , defined semi-circular shell towards the south-east and weak dif- 
fuse emission towards the north-west. Optical observations re- 
vealed strong [OIII] filaments, which generally coincide with CTA 1 was argued to be a site of enhanced interstellar plasma 
the radio shell (Fesen et al. 1981, 1983; Mavromatakis et al. turbulence by analysing the structure functions for rotation mea- 
2000). Strong centre-filled X-ray emission was also observed, sures (RMs) of extragalactic sources in the line-of- sight through 
which consists of both a thermal component heated by the re- and outside the SNR (Simonetti 1992). In principle it would be 
verse shock and a non-thermal component powered by a then straightforward to compare the RM structure function of extra- 
putative pulsar (Seward et al. 1995; Slane et al. 1997). The pul- galactic sources directly with that of the SNR. Polarised emis- 
sar together with its wind nebula was recently detected by Fermi s ion has been observed at 2695 MHz (Sieber et al. 1981) and 
(Abdo et al. 2008). All observations indicate that CTA 1 is a 1720 MHz (Sieber et al. 1979). However, measurements at these 
composite SNR in the Sedov phase. two bands cannot warrant a reliable RM map for the SNR be- 
CTA 1 is a text-book example for the breakout phe- cause of depolarisation at 1720 MHz. Polarisation observations 
nomenon (Pineault et al. 1997), similar to the SNR VRO 42.05.1 at higher frequencies such as 4.8 GHz are therefore essential to 
(Landecker et al. 1982). The blast wave of the supernova sweeps calculate RMs. 
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CTA 1 resides at / = 119°5, b = Mf.2 and experiences little 
obscuration by the emission from the plane at that high latitude. 
This makes it a well- suited object to study the spatial variation 
of the spectral index. In this paper we present new A6 cm and 
All cm observations of CTA I, which is one of the targets of 
the campaign to map large northern sky SNRs at A6 cm. This 
project is part of the Sino-German A6 cm polarisation survey of 
the Galactic plane (Sun et al. 2007; Gao et al. 2010; Sun et al. 
2011; Xiao et al. 2011; Gao et al. 2011). Results on some large 
SNRs have been already reported, such as G65.2+5.7 (Xiao et al. 
2009), the Cygnus Loop (Sun et al. 2006), HB 3 (Shi et al. 2008), 
G156.2+5.7 (Xu et al. 2007), and S 147 (Xiao et al. 2008). 

We describe the new A6 cm and All cm observations of 
CTA 1 as well as other data used, and present the maps in Sect. 2. 
The results and discussions on the spectrum of integrated flux 
densities, maps of spectral index and RM are given in Sect. 3. 
Particularly the detection of a large high-latitude Faraday screen 
is reported. We summarize our results in Sect. 4. 

Throughout the paper we adopt the distance to the SNR of 
1.4 + 0.3 kpc determined from H I measurements (Pineault et al. 
1993). 

2. New observations and other data used 

2.1. Urumqi A6 cm observations 

CTA 1 was observed with the Urumqi 25 -m telescope of 
National Astronomical Observatories, Chinese Academy of 
Sciences, between September and December 2004. Details on 
the A6 cm receiving system were already given by Sun et al. 
(2006, 2007). A field of 3° x 3° size, centred at RA = h 15 m , 
Dec = 72°48' (unless otherwise noted, all the Equatorial coor- 
dinates are in the epoch of J2000.0 hereafter), was observed in 
raster- scans along both Right Ascension and Declination direc- 
tions. All relevant observational parameters are listed in Table 1. 

Table 1. Observation parameters at A6 cm and All cm for 
CTA 1. 



Wavelength 


A6 cm 


All cm 


Frequency 


4.8 GHz 


2.639 GHz 


Bandwidth 


600 MHz 


80 MHz 


Resolution 


9:5 


4:4 


T 

1 sys 


22 K 


17 K 


T B [K]/S [Jy] 


0.164 


2.52 


Scan velocity 


2°/min 


3°/min 


Sub- scan separation 


3' 


2' 


Coverages 


9 


4 


rms of / (07) 


0.5 mK 


3 mK 


rms of PI (cr PI ) 


0.3 mK 


1 mK 


Primary calibrator 


3C 286 


3C 286 


Flux density 


7.5 Jy 


10.4 Jy 


Polarisation percentage 


11% 


9.9% 


Polarisation angle 


33° 


33° 



The data processing procedures were described by e.g. Sun 
et al. (2007) and Gao et al. (2010). The raw data from the obser- 
vations are maps of Stokes /, U, and Q stored in NOD2 format 
(Haslam 1974). For each individual map, spikes were removed, 
the baselines were adjusted, and the scanning effects were sup- 
pressed using the method developed by Sofue & Reich (1979). 
The edited map was further multiplied by a calibration factor de- 
termined from 3C 286 to convert map-units into units of main- 
beam brightness temperature. Maps observed in orthogonal di- 



rections were then added in the Fourier domain (Emerson & 
Grave 1988) to further eliminate scanning effects and yield the 
final results. Instrumental polarisation was cleaned according to 
observations of the unpolarised calibrator 3C 295 following the 
method described by Sun et al. (2007). 

We show the new A6 cm total intensity and polarisation 
maps in Fig. 1. The polarisation intensity (PI) was obtained 

as PI = ^]U 2 + Q 2 -l.2cr 2 U Q following Wardle & Kronberg 

(1974) to correct for the positive noise offset. The polarisation 
angle (iff) was calculated as if/ = |atan^. Here crjj and ctq are 
the rms-noise of U and Q. The rms-noise is 0.5 mK T B for /, 
and 0.3 mK T B for U, Q, and PI, which are quite close to the 
theoretical expectations. 

Similar to previous radio observations of CTA 1 (Sieber 
et al. 1981; Pineault et al. 1997), a very bright semicircular shell 
roughly centred at RA « h 6 m , and Dec « 72°48' with a ra- 
dius of about 50' clearly outlines half of the SNR. The oppo- 
site half of CTA 1 mainly manifests weak and diffuse emis- 
sion, which characterizes the "breakout" region. A bright arc 
protruding from the southern shell is clearly seen, which is 
called central branch hereafter. This corresponds to the "bridge" 
named by Pineault et al. (1997). The weak hooklike feature (just 
above the gamma-ray pulsar indicated in Fig. 1) as reported 
by Pineault et al. (1997) connecting the top parts of the cen- 
tral branch and eastern shell can be identified. The strong com- 
pact source coinciding with the shell is planetary nebula NGC 40 
(RA = h 13 m 2 s , Dec = 72°30'31"), discussed in Appendix A. 

Strong polarised emission corresponding to the shell and the 
central branch was detected. The B -vectors follow the shell and 
the central branch. As shown below in Sect. 3.3.1, the intrinsic 
RM of CTA 1 is up to 40 rad m~ 2 . This means the polarization 
angles at A6 cm deviate less than 9° from their intrinsic values. 
Thus the B -vectors at A6 cm shown in Fig. 1 fairly represent the 
orientation of the magnetic field. Near NGC 40, the polarisa- 
tion is much weaker than in the neighbouring area, which splits 
the shell into two parts (the eastern shell and the southern shell 
hereafter). This low-polarisation region is probably caused by 
a dense cloud pre-dating the supernova explosion, which drasti- 
cally decreases the particle acceleration efficiency (Pineault et al. 
1997). There is polarised emission connecting the central branch 
and the eastern shell. The average polarisation percentage (PC) 
at A6 cm is about 37% for the eastern shell, about 30% for the 
central branch, and about 28% for the southern shell. 

2.2. Effelsberg A11 cm observations 

CTA 1 was observed by Sieber et al. (1979) using the Effelsberg 
100-m telescope at 2695 MHz. The system temperature was 
about 100 K at that time. A new All cm receiver was installed in 
2005, which has a much lower system temperature of 17 K. We 
conducted observations of CTA 1 at All cm using the new sys- 
tem in November and December 2005 and January 2006. Four 
coverages were obtained. The relevant observation parameters 
are listed in Table 1 . 

The data were processed following standard procedures for 
Effelsberg continuum observations, which were the same used 
to reduce the Urumqi A6 cm observations. 

The total intensity and polarisation maps of CTA 1 at Al 1 cm 
are shown in Fig. 2. Both the total intensity and the polarisation 
maps resemble the corresponding A6 cm maps, but have a higher 
angular resolution. The breakout region can be tracked further to 
the north than in the A6 cm map, indicating a very steep spec- 
trum. 
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Fig. 1. Urumqi A6 cm maps. Total intensities are displayed as 
image and contours in the upper panel. Polarisation intensities 
are displayed as image in the lower panel, overlaid with bars 
showing B-vectors (i.e. observed E-vectors rotated by 90°) and 
total intensity contours. The lengths of the bars are proportional 
to polarisation intensities with a cutoff below 5xcr PI . Contours 
start at 5 mK 7b and run in steps of 10 mK T B . The gamma-ray 
pulsar discovered by Fermi is marked as a star. The eastern shell, 
central branch and southern shell are sketched by blue dashed 
lines. 



A number of details can be revealed in the polarisation 
image. The northern end of the eastern shell and the central 
branch seem to split CTA 1 into two branches. Towards the 
low-polarisation region near NGC 40, we note several discrete 
small patches as the one at RA = h ll m , Dec = 72°24', which 
is barely visible in the A6 cm polarisation image. These could 
be fragments of a dense cloud disrupted after the passage of 
the shock wave from the supernova explosion as discussed by 
Pineault et al. (1997). 

The average polarisation percentage is about 30% for the 
eastern shell, 22% for the central branch, and 26% for the south- 
ern shell. 
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Fig. 2. Same as Fig. 1 but for Effelsberg Al 1 cm maps. The start- 
ing level is 30 mK T B and the interval is 80 mK T B for the con- 
tours. 



2.3. 408 MHz and 1420 MHz total intensity data 

408 MHz and 1420 MHz maps were published by Pineault 
et al. (1997), which were observed with the DRAO synthesis 
array, where the missing large-scale components were included 
from the 408 MHz all-sky survey (Haslam et al. 1982) and 
1420 MHz observations made with the Effelsberg telescope in- 
cluding data from the 1420 MHz Stockert survey (Reich 1982) 
for the broadest structures. The maps have angular resolutions of 
3f5 at 408 MHz and V at 1420 MHz. We use these data to con- 
struct a spectral index map together with the A6 cm and Al 1 cm 
observations. 



2.4. Effelsberg 1400 MHz polarisation data 

We extracted 1400 MHz polarisation data for the area of CTA 1 
from an unpublished section of the Effelsberg Medium Latitude 
Survey (EMLS), which was described by Uyaniker et al. (1998) 
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and Reich et al. (2004). The resolution is 935. The 1420 MHz 
total intensity map from the combined DRAO and Effelsberg ob- 
servations is of high quality (Pineault et al. 1997) and therefore 
not replaced by total intensity data from the EMLS. 

The A21 cm polarisation map from the EMLS is shown in 
Fig. 3. The polarisation distribution along the eastern shell is 
quite similar to that observed at A6 cm and All cm. For the 
southern shell, eastwards of RA « h 4 m , however, there is al- 
most complete depolarisation. Parts of the central branch also 
experience large depolarisation. The average polarisation per- 
centage is about 26% for the eastern shell, 14% for the cen- 
tral branch, and 10% for the southern shell. Polarisation patches 
originating within the interstellar medium surround CTA 1 and 
most likely also exist along the line-of-sight to the SNR. They 
cause an overestimate of the A21 cm polarisation percentage of 
CTA 1. 




Right Ascension (J2000) 

Fig. 3. Same as Fig. 1 (lower panel) but for the EMLS A21 cm 
map. Total intensity contour levels start at 200 mK 7b and run 
in steps of 200 mK T^. 



3. Spectrum and rotation measure of CTA 1 

3.1. Integrated flux densities and spectrum of CTA 1 

We used flux-density measurements of CTA 1 at 408 MHz, 
1420 MHz, 2639 MHz, and 4800 MHz to constrain the spec- 
trum of this SNR. Data at lower frequencies compiled by Sieber 
et al. (1981) were not used, because their very low angular res- 
olution makes it difficult to assess the contribution of compact 
sources. 

Compact, point like sources visible in the total intensity 
maps at 408 MHz, 1420 MHz, and 2639 MHz, were fitted by 
a 2D-elliptical Gaussian to subtract them. Most of the extra- 
galactic sources are very weak at A6 cm and confused with 
CTA 1 . Therefore their flux densities cannot be measured from 
the A6 cm map, and extrapolation from low frequency observa- 
tions is the only way to estimate their contribution. We obtained 
a map of sources at 1420 MHz and scaled it to 4.8 GHz using a 
spectral index of a = -0.9 (e.g. Zhang et al. 2003). The spectral 
index is defined as S v oc y a with S v being the flux density at 



frequency v. The scaled map was then smoothed to 9'. 5 and sub- 
tracted from the A6 cm map. NGC 40 has been treated separately 
according to its spectrum shown in Fig. A.l. A hyper-plane cal- 
culated from average values from the four map-corners was then 
subtracted from each map to set the surroundings of the SNR to 
zero. 

We measured the integrated flux density of CTA 1 to be 
11.6+1.2 Jy at A6 cm and 20.3+2.0 at All cm after discounting 
the contribution of extragalactic sources. We corrected the flux 
densities obtained by Pineault et al. (1997) accordingly and used 
integrated flux densities of 31+3 Jy and 60+4 Jy at 1420 MHz 
and 408 MHz, respectively. Fitting the flux density values at 
these four frequencies yields a spectral index of a = -0.63+0.05 
(Fig. 4), consistent with that reported by Pineault et al. (1997) of 
a = -0.57 + 0.006. 



CD 
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Fig. 4. Fitted spectrum of integrated flux densities. 



As a crosscheck for the integrated flux density spectral index, 
we made TT-plots (Turtle et al. 1962) between A6 cm and the 
other three frequencies (Fig. 5). With this method the influence 
of a possible incorrect zero-level setting can be controlled. The 
TT-plot spectral index (J3) from brightness temperatures relates 
to a as a = ft + 2. As shown in Fig. 5, the spectral index from 
TT-plots is f3 = -2.58 + 0.03 between A6 cm and 408 MHz, = 
-2.61 + 0.03 between A6 cm and A21 cm, and p = -2.61 + 0.05 
between A6 cm and All cm. These values convincingly confirm 
the spectral index we obtained from the integrated flux densities. 
The group of outliers at 40-50 mK 7b in the TT-plots in Fig. 5 
are from the residuals when subtracting NGC 40 from the A6 cm 
map. 



3.2. Spectrum steepening 

A steepening of the spectrum from the shell towards the break- 
out region was already reported by Pineault et al. (1997) based 
on the 408 MHz and 1420 MHz observations. The new A6 cm 
and All cm measurements allow us to study the spatial varia- 
tions of the spectral index more accurately, because of the wider 
frequency range. 

Towards the SNR, the observed intensity is a sum of the in- 
trinsic emission from the SNR and the fluctuations of the diffuse 
interstellar medium on degree scales. To set up a common zero- 
level for all the maps, an offset has to be applied for each fre- 
quency before calculating a spectral index map of CTA 1 . The 
source-removed maps at 408 MHz, 1420 MHz and 2639 MHz 
were smoothed to an angular resolution of 9f5. We assume the 
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Fig. 7. Spectral index map of CTA 1 calculated from the four 
maps between 408 MHz and 4.8 GHz. The contours indicate the 
smoothed Al 1 cm total intensities shown in Fig. 6. 



offset to be zero at 4.8 GHz. This is reasonable because the large- 
scale emission is very weak at the high-latitude of CTA 1 . We 
then added or subtracted a constant offset from the maps at each 
of the other three frequencies to assure the TT-plots versus A6 cm 
intersect at zero brightness temperature (Fig. 5). Total intensities 
of 830 mK, 87 mK and 24 mK were subtracted from maps at 
408 MHz, 1420 MHz and 2639 MHz, respectively, to properly 
settle the zero-levels. The adjusted maps are displayed in Fig. 6. 

For each pixel of the region encompassing CTA 1 from zero- 
level corrected maps at the four frequencies (Fig. 6), we made a 
linear fit for intensities versus frequencies in logarithmic scale to 
obtain the spectral index map. For regions at Dec > 73°30 / , only 
intensities at 408 MHz, 1420 MHz and 2639 MHz were used 
since no diffuse emission from the SNR could be detected at 
4.8 GHz exceeding the noise level. The result is shown in Fig. 7. 
The errors are generally less than 0.1 and become larger towards 
the very northern part of CTA 1, where the spectral index de- 
creases to about a = -I. 

Spectral index values are between a = -0.5 and a = -0.65 
towards the shell and the central branch, and gradually become 
smaller towards the breakout region. The variation is up to about 
Aa = 0.3, which confirms the results by Pineault et al. (1997). 

The mechanism for the spectrum steepening is not yet clear. 
It is widely accepted that cosmic-ray electrons experience diffu- 
sive shock acceleration (DSA) in SNR shock fronts and become 
relativistic. The spectral index of the synchrotron emission ob- 
served from these electrons can be written as a = -3/2(r - 1), 
where r is the shock compression ratio. To produce the observed 
spectral index, the compression ratio is about 3.5-4 for the shell, 
and 2.7-2.9 for the breakout region. A decrease of the com- 
pression ratio requires a decreasing shock Mach number, which 
is possible in case of a low gas density in the breakout region 
(Pineault et al. 1997). However, whether DSA is still important 
in the breakout region remains very uncertain. 

It could be the case that the steep- spectrum emission of the 
breakout region originates from higher-energy electrons than in 
the SNR shell because of a weaker magnetic field there. These 
high-energy electrons may have a steeper spectrum than the low- 
energy electrons due to synchrotron aging. They may have been 



accelerated and diffused out from the SNR shell. In this case 
a steepening at higher frequencies should be seen in the SNR 
spectra. However, no indication of a spectral break is found in 
the frequency range between 408 MHz and 4800 MHz, which 
means that any spectral break due to synchrotron aging should 
be at higher frequencies. Sensitive observations of CTA 1 at even 
higher frequencies are needed including measurements of the ex- 
tended diffuse emission from the breakout region, although it is 
difficult to map such a large object with arcmin angular resolu- 
tion. 

3.3. RM map of CTA 1 
3.3.1. RM determination 

We smoothed the Al 1 cm U and Q data to an angular resolution 
of 9' 5 and derived a map of polarisation angles. According to 
the polarisation angles at A6 cm and All cm, we calculated RMs 

as, 

^6 cm - <Al 1 cm 



RM 



6 cm 



" ^11 cm 



(1) 



RM is proportional to the integral of the thermal electron den- 
sity multiplied by the magnetic field parallel to the line-of- sight. 
Pixels with brightness temperatures less than 5 x cr PI were not 
included. The results are shown in Fig. 8. RMs of extragalac- 
tic sources in the CTA 1 area were taken from the catalogue by 
Taylor et al. (2009) and displayed in Fig 8. Their RMs range be- 
tween ±40 rad m~ 2 and are much smaller than the RM ambiguity 
of 348 rad m~ 2 between A6 cm and All cm. Therefore we always 
use the minimum polarisation angle difference to calculate RM. 

A 20 rad m~ 2 



74° 



73° 



72° 




n 20 ] 



n 12 ! 



n 04= 



23 h 56 m 



Right Ascension (J2000) 



Fig. 8. RM contours overlaid on the A6 cm polarisation image. 
The contour levels are (red), 10 (magenta), 20 (orange), and 
30 rad m~ 2 (red+magenta) for positive RMs and -10 (blue), 
-20 (cyan), -30 (green+yellow), and -40 rad m~ 2 (blue+cyan) 
for negative RMs. The filled (open) circles represent positive 
(negative) RMs of extragalactic sources taken from Taylor et al. 
(2009). The dashed line marks Galactic latitude b = 9° and 
b = 11°. 

The RM distribution of CTA 1 exhibits a clear pattern. 
Towards the eastern shell and the central branch, RMs are al- 




Fig. 5. TT-plots between A6 cm (4800 MHz) and other frequencies as indicated in the panels. 





Right Ascension (J2000) Right Ascension (J2000) 

Fig. 6. Total intensity maps at four frequencies. The sources have been subtracted and the zero-levels have been corrected. All the 
maps have a common resolution of 9f5. 



ways positive, and always negative along the southern shell. The 
average of RMs for the eastern shell and the central branch is 
about +10 rad m~ 2 and +15 rad m~ 2 . For the southern shell, the 
absolute values of RMs are as large as about 50 rad m~ 2 . The 
RM map by Sieber et al. (1981) based on the Effelsberg Al 1 cm 
and A 18 cm observations qualitatively agrees with our result. 
Interestingly the RMs of extragalactic sources show a similar 
pattern as CTA 1 , but extend to a much larger region as discussed 
in Sect. 3.3.2. 



3.3.2. A Faraday screen in front of CTA 1 

The RM pattern of CTA 1 as well as that of extragalactic sources 
(Fig. 8) strongly suggests, that there is a Faraday screen in the 
direction of the eastern shell and the central branch of CTA 1. 
The Faraday screen exceeds the size of CTA 1 as it can be 
clearly seen from a larger field containing RMs from extragalac- 
tic sources (Fig. 9). Therefore the Faraday screen should be lo- 
cated in front of CTA 1, which has a distance of about 1.4 kpc. 
Its centre is roughly at / = 121°,/? = 11° with a poorly con- 
strained size of about 3°. To estimate the RM of the Faraday 
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screen, we first need to investigate RMs from the SNR and the 
Galactic foreground. 



• 20 rad m 
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Galactic Longitude 
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Fig. 9. Similar to Fig. 8 but for a larger field in Galactic coordi- 
nates. The contours outline positive (+10 rad m~ 2 , red) and neg- 
ative (-10 rad m~ 2 , blue) RMs. The circle marks the approximate 
boundary of the Faraday screen. 

The intrinsic RMs of CTA 1 can be estimated from depo- 
larisation. Polarised emission originating from different depth 
experiences different Faraday rotation, and summing up all the 
emission components along the line-of-sight may partly cancel 
polarised emission components. We define the relative depolari- 
sation (DP) at Al 1 cm and A21 cm as DP = PC/PC A6 cm , where 
PC is the polarisation percentage. The All cm and A21 cm de- 
polarisation maps are shown in Fig. 10. Following Sokoloff et al. 
(1998) wavelength-dependent depolarisation is related to RM as, 



Dp = sin(2|RM|^ 2 ) x 4 
A sin(2|RM|^ 2 ) A 2 ' 



(2) 



where we take Aq = 6.25 cm (4.8 GHz). 

The average depolarisation is about 0.9, 0.8, and 0.8 for the 
eastern shell, the central branch, and the southern shell, respec- 
tively. This implies a small amount of depolarisation at Al 1 cm. 
At All cm, the average depolarisation is about 0.9 for the eastern 
shell, which means an intrinsic absolute RM of 10 rad m~ 2 . The 
average depolarisation of about 0.6 towards the central branch 
requires an absolute RM value of of about 18 rad m~ 2 . Towards 
the southern shell except for several fragments, there is complete 
depolarisation at A21 cm. The average depolarisation of the frag- 
ments is about 0.3, implying an absolute RM value larger than 
about 25 rad m~ 2 . 

The intrinsic RMs of CTA 1 are negative following the reg- 
ular Galactic magnetic field direction in this area. The observed 
positive RMs towards the eastern shell and the central branch are 
caused by the Faraday screen in front of CTA 1 . 

Based on the 3D-emission models of the Milky Way by Sun 
et al. (2008) and Sun & Reich (2010), a RM from the diffuse 
interstellar medium up to the 1 .4 kpc distance of CTA 1 was cal- 
culated to be about -17 + 8 rad m~ 2 . Unfortunately the newly dis- 
covered gamma-ray pulsar by Fermi (Abdo et al. 2008) does not 
yet have a radio detection and thus a measured RM. However, 
the pulsar B0105+68, about 6° offset from the centre of CTA 1 
has a distance of about 2.6 kpc and a RM of -46 rad m~ 2 (Mitra 
et al. 2003). A linear interpolation to 1.4 kpc distance yields a 
RM of -24 rad m~ 2 , consistent with that from the 3D-emission 
models. 
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Fig. 10. Relative depolarisation relative to A6 cm at Al 1 cm (up- 
per panel) and A21 cm (lower panel). 



The RM of the Faraday screen can be calculated from the ob- 
served RMs towards the eastern shell and the central branch of 
CTA 1 by subtracting the intrinsic CTA 1 values and the Galactic 
foreground contribution. The result is about 37-50 rad m~ 2 . In 
the following we use a mean value of about 45 rad m~ 2 for 
the Faraday screen. For the southern shell, which is not influ- 
enced by the Faraday screen, the intrinsic RMs plus that from 
the Galactic foreground roughly agree with the observed ones. 

The Faraday screen will cause a rotation of the polarisation 
angles at A21 cm by more than 100°. However, we cannot find 
a correspondence with the EMLS data. This implies that most 
of the polarised emission observed at A21 cm originates from 
regions closer than the Faraday screen. 

Assuming the Faraday screen to be spherical with a mean 
size of about 3°, the regular magnetic field para llel to the line- 
of-sight can be estimated asB\\ = 3.6 RM/ ^I\\ a d. Here Iu a is the 
Ha intensity in R, d is the distance in pc, the electron tempera- 
ture is assumed to be 8000 K, and the extinction is neglected. 
There is no enhanced Ha emission visible towards the Faraday 
screen in the Wisconsin H- Alpha Mapper northern sky survey 
(Haffner et al. 2003). The average and the fluctuations of the 
Ha intensity are about 4 R and 1 R, respectively, which implies 
an upper limit of about 7 R for the Faraday screen. For an as- 
sumed distance of 500(1000) pc, the depth of the screen along 
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line-of-sight is about 26(52) pc and the lower limit of the regular 
magnetic field parallel to the line-of-sight is about 2.7(1.9) fiG. 

3.3.3. RM structure functions 

Simonetti (1992) found that CTA 1 is a site of enhanced interstel- 
lar plasma turbulence based on a structure function analysis for a 
number of RMs of extragalactic sources at line-of-sights through 
and outside the SNR. The RM map we obtained above allows us 
to compare the structure function of RMs from CTA 1 directly 
with that of RMs from extragalactic sources by Simonetti (1992) 
and Taylor et al. (2009). 

The second order structure function is a frequently used tool 
to infer turbulent properties of the interstellar medium (e.g. Sun 
& Reich 2009; Stil et al. 201 1), which is defined as, 

D(S6) =< [RM(0) - RM(0 + 80)f > . (3) 

Here 80 is the angular separation and < • • • > stands for the en- 
semble average. We calculated the structure functions from the 
RM map of CTA 1 for the southern shell, where RM is negative, 
and for the eastern shell and the central branch, where RM is 
positive (see Figs. 8 and 9). We also obtained the structure func- 
tion for RMs of extragalactic sources shown in Fig. 9. Note that 
several sources coinciding with CTA 1 (Fig. 9) were excluded in 
order to compare with the RM structure function for extragalac- 
tic sources with line-of-sights outside of CTA 1 by Simonetti 
(1992). The results are shown in Fig. 11. 
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Fig. 11. Structure functions for RMs of extragalactic sources 
taken from Taylor et al. (2009), the southern shell of CTA 1 (neg- 
ative RMs), and the eastern shell and the brigde of CTA 1 (posi- 
tive RMs). The squared RM differences for extragalactic sources 
at line-of-sights through and outside CTA 1 from Simonetti 
(1992) are also shown. 

The structure functions for both positive and negative RMs 
of CTA 1 can be well fitted by a power law with a slope of about 
0.85 in the range of 10' < 80 < 30', where 10' reflects the angu- 
lar resolution and 30' is about the width of the shell. Note that 
the RMs of CTA 1 consist of contributions from the SNR and the 
Galactic foreground for the southern shell. For the eastern shell 
and the central branch, the contribution from the Faraday screen 
adds, which does not change the slope of the structure function. 
This indicates that the magnetic field in the Faraday screen along 
the line-of-sight is very regular, so that no additional fluctuations 
were introduced. The different amplitudes of the two structure 
functions stem from the different absolute values of RM, which 
are also reflected in the depolarisation properties (see Sect. 4.2). 



The squared RM differences of extragalactic sources at line- 
of-sights through and outside CTA 1 from Simonetti (1992) are 
also shown in Fig. 11. At larger scales of 80 > 20', the ampli- 
tudes of RM structure functions for all extragalactic sources are 
similar. At smaller scales of 80 < 3', the squared RM differ- 
ences for extragalactic sources at line-of-sights through CTA 1 
are larger than outside CTA 1. Therefore Simonetti (1992) pro- 
posed that CTA 1 induces enhanced plasma turbulence with an 
outer scale in the range of 3 / -20 / . As indicated in Fig. 11, the 
amplitude of the RM structure function for the shell and central 
branch regions of CTA 1 roughly agrees with that for extragalac- 
tic sources at line-of-sights outside CTA 1, when extrapolated to 
smaller scales. This implies that the scenario of enhanced plasma 
fluctuation proposed by Simonetti (1992) is limited at scales be- 
low about 3'. Clearly more data are needed to settle this issue. 

4. Summary 

We performed new A6 cm and All cm continuum observations 
of the shell-type SNR CTA 1 . We obtained integrated flux densi- 
ties of 11.6+1 .2 Jy at A6 cm and 20.3+2.0 Jy at Al 1 cm after sub- 
traction of the contributions from compact extragalactic sources. 
Together with values at 408 MHz and 1420 MHz we derived a 
spectral index of a = -0.63 + 0.05, which was confirmed by 
TT-plots. We also obtained a spectral index map based on data 
at 408 MHz, 1420 MHz, 2639 MHz, and 4800 MHz. The map 
shows a spectral index variation up to about Aa = 0.3 from the 
shell towards the breakout region. The reason for the spectrum 
steepening needs further investigation. 

We obtained a RM map using polarisation data at A6 cm and 
All cm. A clear pattern of negative RMs towards the southern 
part and positive ones towards the eastern part can be seen from 
the RM map. We argued that there exists a Faraday screen of 
roughly 3° in size located in front of CTA 1 covering its north- 
eastern part. The RM of the Faraday screen is about 45 rad m~ 2 . 
Assuming a distance of 500(1000) pc for the Faraday screen the 
lower limit for the regular magnetic field parallel to the line-of- 
sights is 2.7(1.9) yuG. 

We calculated RM structure functions for different parts of 
CTA 1 and could not find an influence by the foreground Faraday 
screen, which indicates a very regular enhanced magnetic field. 
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Fig. A.l. Spectrum of integrated flux densities of NGC 40. The 
integrated flux densities collected by Sieber et al. (1981) are 
shown by squares. More flux densities were reported at 330 MHz 
(Rengelink et al. 1997), 1420 MHz (Condon et al. 1998; Pineault 
et al. 1993), 4.85 GHz (Becker et al. 1991), 30 GHz (Pazderska 
et al. 2009), and 43 GHz (Umana et al. 2008) and shown by open 
circles. The new A6 cm and All cm measurements are marked 
as filled circles. 



in Fig. A.l. The two-component radio emission model of plane- 
tary nebulae by Siodmiak & Tylenda (2001) was used to fit the 
spectrum, which yields an opacity of r = 0.0054(v/5 GHz) -21 . 



Appendix A: NGC 40 

NGC 40 (RA = h 13 m 2 s , Dec = 72°30'31") is a well-known 
planetary nebula where its central star shows the Wolf-Rayet 
phenomenon (Grosdidier et al. 2001). Its distance is very uncer- 
tain and estimates range from 0.5 kpc to 3.5 kpc (Phillips 2005). 
From the polarisation images at A6 cm (Fig. 1), Al 1 cm (Fig. 2), 
and A21 cm (Fig. 3), it is evident that NGC 40 does not have 
an effect on polarisation distribution and cannot cause the low- 
polarisation region, which divides CTA 1 into an eastern and a 
southern shell. Therefore its distance should be larger than that 
of CTA 1 unless it causes no Faraday rotation. 

The integrated flux density of NGC 40 was measured to be 
595 ± 47 mJy at A6 cm and 583 ± 58 mJy at Al 1 cm. Besides the 
flux density measurements collected by Sieber et al. (1981) we 
found some new ones from literature and plotted both datasets 



